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1  INTRODUCTION 

There  is  currently  a  great  deal  of  interest  in  the  use  of  polarimetry  for  radar  remote  sensing.  In  this  context, 
different  and  important  objectives  are  to  classify  Earth  terrain  components  within  a  fully  polarimetric  SAR 
image  and  then  extract  physical  information  from  the  observed  scattering  of  microwaves  by  surface  and 
volume  structures.  The  most  important  observable  measured  by  such  radar  systems  is  the  3x3-coherency 
matrix  [J].  This  matrix  accounts  for  local  variations  in  the  scattering  matrix  and  is  the  lowest  order  operator 
suitable  to  extract  polarimetric  parameters  for  distributed  scatterers  in  the  presence  of  additive  (system) 
and/or  multiplicative  (speckle)  noise.  In  the  first  part  of  this  paper,  the  most  important  Target  Polarimetry 
descriptors:  Sinclair  Matrix,  target  vectors,  coherency  matrix  and  the  covariance  matrix  as  well  are 
presented,  their  interconnections  and  equivalences  will  be  shown  together  with  the  respective 
transformations. 

Speckle  appearing  in  synthetic  aperture  radar  (SAR)  images  is  due  to  the  coherent  interference  of  waves 
reflected  from  many  elementary  scatterers  and  causes  degradation  and  makes  automatic  image  segmentation 
and  scene  description  difficult.  The  speckle  reduction  problem  is  more  complicated  for  polarimetric  SAR 
than  a  single  polarization  SAR,  because  of  the  difficulties  of  preserving  polarimetric  properties  and  of 
dealing  with  the  cross-product  terms.  The  first  part  of  this  paper  is  ended  by  a  presentation  and  a  description 
of  polarimetric  speckle  filters  preserving  polarimetric  properties  and  statistical  correlation  between  channels, 
not  introducing  crosstalk,  and  not  degrading  the  image  quality.  The  impact  of  using  this  polarimetric  speckle 
filtering  on  terrain  classification  is  quite  dramatic  in  boosting  classification  performance. 

Many  targets  of  interest  in  radar  remote  sensing  require  a  multivariate  statistical  description  due  to  the 
combination  of  coherent  speckle  noise  and  random  vector  scattering  effects  from  surface  and  volume.  For 
such  targets,  it  is  of  interest  to  generate  the  concept  of  an  average  or  dominant  scattering  mechanism  for  the 
purposes  of  classification  or  inversion  of  scattering  data.  Target  Decomposition  theorems  are  aimed  at 
providing  such  an  interpretation  based  on  sensible  physical  constraints  such  as  the  average  target  being 
invariant  to  changes  in  wave  polarization  basis.  Among  the  existing  Polarimetric  Target  Decomposition 
theorems  -  coherent  (Krogager,  Cameron  ...),  non-coherent  (Huynen,  Barnes  ...)  model-based  decomposition 
(Freeman)  or  eigenvector-based  decomposition  (Cloude,  Van  Zyl)  -  the  H/A/a  Decomposition  Theorem, 
proposed  by  S.R.  Cloude  and  E.  Pottier  in  1997  for  extracting  average  parameters  from  experimental  data  is 
presented  and  discussed.  Based  on  an  eigenvalues  analysis  of  the  coherency  matrix,  this  decomposition 
theorem  employs  a  3 -level  Bernoulli  statistical  model  to  generate  estimates  of  the  average  target  scattering 
matrix  parameters. 


2  POLARIMETRIC  TARGET  DESCRIPTORS 


2.1  Target  vector  for  backscatter  problems. 

Radar  Polarimetry  (Polar:  polarization  Metry:  measure)  is  the  science  of  acquiring,  processing  and  analysing 
the  polarization  state  of  an  electromagnetic  field.  The  polarization  information  contained  in  the  waves 
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backscattered  from  a  given  medium  is  highly  related  to  its  geometrical  structure  and  orientation  as  well  as  to 

its  geophysical  properties  such  as  humidity,  roughness  and  conductivity  of  soils, . 

Radar  Polarimetry  deals  with  the  full  vector  nature  of  polarized  electromagnetic  waves,  and  when  the  wave 
passes  through  a  medium  of  changing  index  of  refraction,  or  when  it  strikes  an  object  or  a  scattering  surface 
and  it  is  reflected;  then,  characteristic  information  about  the  reflectivity,  shape  and  orientation  of  the 
reflecting  body  can  be  obtained  from  the  2x2  coherent  backscattering  matrix  or  Sinclair  matrix  [Boemer 
1998]. 

An  important  development  in  our  understanding  of  how  to  best  extract  physical  information  from  the 
classical  2x2  coherent  backscattering  matrix  [5]  has  been  achieved  through  the  construction  of  system 
vectors  [Cloude  1996][Cloude  1997].  We  represent  this  vectorization  of  a  matrix  by  the  vector  V(.)  built  as 
follows  : 


[5] 


k  =  V{[S])  =  -2Tracei^S][¥^ 


(1) 


where  Trace([A])  is  the  sum  of  the  diagonal  elements  of  matrix  [A]  and  [  y/\  is  a  set  of  2x2  complex  basis 
matrices  which  are  constructed  as  an  orthonormal  set  under  an  hermitian  inner  product  [Cloude  1986][ 
Cloude  1996].  There  exist  in  the  literature  different  basis  sets,  but  the  special  set  used  to  generate  3x3 
coherency  matrix  [7]  is  based  on  linear  combinations  arising  from  the  Pauli  matrices,  and  is  given  by: 


42 


1  0 
0  1 


(2) 


where  the  factor  arises  from  the  requirement  to  keep  the  norm  of  the  target  vector  k  invariant,  equal  to 
the  Frobenius  norm  (Span)  of  the  backscattering  matrix  [*S],  namely  the  total  power  scattered  by  the  target. 
The  target  vector  k  has  the  explicit  form  shown  in  (3). 


k  = 


(3) 


With  such  a  vectorization  we  can  then  generate  a  coherency  matrix  from  the  outer  product  of  the 
target  vector  k  with  its  conjugate  transpose.  For  the  monostatic  case,  the  3x3  hermitian  coherency  matrix  [7] 
has  the  following  parameterisation  [Cloude  1996]: 


[T]  =  k-k*T 


2A0  C-jD  H  +  jG 
C  +  jD  B  q  +  B  E  +  jF 
H-jG  E-jF  Bn-B 


(4) 


where  A0,  B0,  B,  C,  D,  E,  E,  G  and  H  are  the  «  Huynen  parameters  »  identified  as  phenomenological  target 
parameters  interrelated  in  such  a  way  as  to  reflect  directly  the  physical  source  of  correlation  effects  in  the 
target.  These  nine  parameters  are  useful  for  general  target  analysis  without  reference  to  any  model,  and  each 
of  them  contains  real  physical  target  information  [Huynen  1970][Pottier  1992]. 

The  parameters  A0,  B0+B  and  B0  -  B ,  called  the  « target  generators  »,  are  connected  respectively  with  the 
symmetry,  irregularity  /  double  bounce  and  non-symmetry  physical  properties  in  the  case  of  a  pure  target 
(strong  scatterer)  [Huynen  1970][Pottier  1992],  or  related  to  surface  scattering,  double-bounce  scattering  and 
volume  scattering  in  the  case  of  a  distributed  target  (natural  media). 

It  is  thus  possible  to  use  these  «  target  generators  »  to  create  a  color  coding  for  PolSAR  images,  by  assigning 
respectively  the  color  red  for  B0  +  B  =  \SHH  -Svv\ ,  green  for  B0-B=\  SHV  | ,  and  blue  for  2A0  =  | SHH  +SW\. 
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For  illustration,  the  well  known  NASA/JPL  AIRSAR  L-Band  full  polarimetric  SAR  image  of  San  Francisco 
Bay  (1988)  and  the  DLR  E-SAR  L-Band  full  polarimetric  SAR  image  of  Oberpfaffenhofen  (Germany)  are 
shown  respectively  on  Figs.  1  and  2. 


Fig.  2:  Optical  image  and  Pauli  color  coded  image  over  Oberpfaffenhofen  (Germany) 


In  these  two  displays,  the  polarimetric  channel  combinations  play  an  important  role  in  separating  ocean 
surface  or  agricultural  areas  (surface  scattering),  forests,  vegetation  and  trees  (volume  scattering),  and  urban 
areas  (double  bounce  scattering)  which  composed  the  different  scenes.  Called  the  “Pauli  coding 
representation”,  this  adopted  color  coding  has  become  today  the  standard  for  PolSAR  image  display. 


There  exists  in  the  literature  a  second  basis  sets,  used  to  generate  the  3x3  covariance  matrix  [C]  and  is  based 
on  a  lexicographic  combinations  of  the  elements  of  the  Sinclair  matrix  and  is  given  by  [Cloude  1996][Cloude 
1997]: 


r 

~1  o' 

,2^2 

'o  r 

'0  0l| 

2 

,  2 

1 

0  0 

0  0 

-°  7Jj 

(7) 


where  the  factor  arises  from  the  requirement  to  keep  the  norm  of  the  target  vector  Q  invariant,  equal  to 
the  Frobenius  norm  (Span)  of  the  backscattering  matrix  [*S],  namely  the  total  power  scattered  by  the  target. 
The  target  vector  /2has  the  explicit  form: 
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(2  = 


S„„  42  S, 


HV  SVv 


1 T 


(8) 


With  such  a  vectorization  we  can  then  generate  a  covariance  matrix  from  the  outer  product  of  the  target 
vector  /2with  its  conjugate  transpose.  For  the  monostatic  case,  the  3x3  hermitian  covariance  matrix  [C]  has 
the  following  parameterisation  [Cloude  1996][Cloude  1997]: 


[c]  =  n-n*T  = 


is 


HH\ 


42  s„x» 

O  O* 

^VV^HH 


42 

2\SHr\‘ 

42  S„S'HV 


C  V 

°HH°W 

42  sm/s; 


Is 


HV^VV 
\2 
w 


(9) 


which  can  be  also  expressed  in  a  well-known  form  given  by  : 


with 


[C]  =  , 


<7=  S 


P  = 


JHH\ 


sHHsvr 


i2wis, 12 


1  p4d  pjy 

P*4s  8  s-^8 

p  Vf  £*  4r$  r 


5=2 


IS 


P  = 


c  c 


Is 


r  = 


Is 


Is 
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^HV^VV 


(10) 


(11) 


The  two  target  vectors  k  and  Q  are  linked  together  with  a  unitary  transformation  which  is  given  by  : 


k  =  -r = 

42 


+1  0  +1 

+1  0  -1 

0  42  0 


Q  <=>  k  =  [u3]n 


(12) 


where  [U3\  is  a  special  unitary  SU(3)  matrix.  It  follows  the  relation  between  the  coherency  matrix  [7]  and  the 
covariance  matrix  [C]: 


[r]  =  [vJ][c][c/3]-'  d3) 

2.2  Elliptical  Polarisation  Basis  Changement  Operation 

In  the  case  of  the  Sinclair  matrix  or  back-scattering  matrix,  the  elliptical  polarisation  basis  changement 
operator  is  given  by  the  following  unitary  congruence  (unitary  con-similarity)  transformation  [Boemer 
1998]: 


[S’]=[U2(t,T,v)]T[Siu4t,T,v)] 


(14) 


where  the  matrix  [ U2 ]  is  given  by  : 
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[u{0,T,v)\  =  [U(<t>)lU(T)][u{vj\ 


cos(^)  -sin(^) 

cos(r)  j  sin(r) 

-  7  V  A 

e  J  0 

sin(^)  cos(^) 

j  sin(  r)  cos(r) 

r\  +  JV 

0  e  J 

where  the  v,  (j),  r  parameters  correspond  to  the  three  geometric  parameters  of  the  polarisation  ellipse 
described  by  the  first  or  principal  Jones  vector  of  the  new  basis. 

The  matrix  [U2\  belongs  to  the  SU(2)  -  Special  Unitary  2x2  matrix  group  constructed  from  the  classical 
unitary  matrices  Pauli  group  given  by  : 


1  0 

cr0  = 

0  |_o  1 

where  the  matrices  verify  err1  =  erf*  and  |det(cr. )|=1 . 

The  group  of  the  special  unitary  matrices  (SU(2))  is  defined  by  [Huynen  1970]: 
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0  -1 
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r  I  J  OCCJ 

[A\=e  p  =cos(a)<r0  +  jsm(a)(T 

It  follows  : 

cos(^)  -sin(^) 
sin(^)  cos(^) 


=cos(^)cr0  -  /sin(^)er3  =e 


-jfa  3 


(16) 


(17) 

(18) 


cos(r) 

./'sin(r) 


jsin(r) 

cos(r) 


=cos(r)cr0  +  ysin(r)cr2 


(19) 


+  jy 

e  ■ 

0 


=cos(  v  )<7n  +  jsin(  v  )<7j 


:+j™, 


These  matrices  verify  the  following  properties: 

[a]  1  =[a]T  and  det[/l]=+  1 


e  J</>(T.i  =e+j0&3 

e+JT(72T  =e+JT(J  2 

+  iv(j1T  +iv<j 
e  J  1  —e 

* 

* 

£+./r<r.  =e~Jxci 

e+Jyai  =e~Jva 

+j(a+fi)ap  _e+ja<yp 

+e+iPa’  ,  „ 

with: 

® p  5  ^  {°"l  5  <72  > 

*pe+J0UT<  =e~JCUT' *p 

(20) 


(21) 

(22) 

(23) 


(24) 


This  special  unitary  elliptical  basis  changement  matrix  can  also  be  described  using  the  parameters  p  and  2 
which  correspond  to  the  polarisation  ratio  of  this  first  or  principal  Jones  vector  of  the  new  basis,  and  are 
given  by  : 
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tan(d>)+  itanvc )  „  (  ( ,\  (  \\ 

p  = -  ,  . — -r\  g  =  v  -  arctanytanyp )tan\z ) ) 

1  -  j  tan\<f>)tan\T) 


(25) 


The  matrix  [ U2 ]  is  thus  given  by  : 


[u,  (S,  r,  v)]  =  [u,  (a  4 )]  =  [u2  ip)\u2  (£ )]  =  -pi 


■(l+jpf 


1  * 
1  -p 

p  1 


ry*  0 


0  e 


+7# 


(26) 


In  the  case  of  the  3x3  coherency  matrix  [7],  the  elliptical  polarisation  basis  changement  operator  is  given  by  : 

(27) 

where  the  matrix  [C/jr]  is  given  by  : 

[U3T  (20,2  t,2v)]  =[U3T  (2  0)\U3T  (2t)][U3T  (2  v)] 


'1 

0 

0 

cos(2r) 

0 

ysin(2r) 

cos(2v) 

-y  sin(2v) 

O' 

0 

cos(2^) 

sin(2^) 

0 

1 

0 

-y  sin(2v) 

cos(2v) 

0 

0 

-sin(2^) 

cos(2^) 

_y'sin(2r) 

0 

cos(2r) 

0 

0 

1 

where  the  v;  (j),  r  parameters  correspond  to  the  three  geometric  parameters  of  the  polarisation  ellipse 
described  by  the  /irsf  or  principal  Jones  vector  of  the  new  basis. 

The  matrix  \U3T\  belongs  to  the  SU(3)  -  Special  Unitary  3x3  matrix  group.  Unfortunately  it  does  not  exist 
any  direct  mathematical  link  between  the  SU(2)  group  matrices  and  the  SU(3)  group  matrices.  To  derive  the 
SU(3)  group  matrices,  we  have  to  deal  with  the  elliptical  polarisation  basis  changement  operation,  and 
identify  for  each  SU(2)  matrix,  its  equivalent  into  the  SU(3)  group.  After  some  derivations,  it  follows: 


[u2(0)\ 


cos(^) 

sin(^) 


-sin(^) 

cos(^) 


>[U3T(20)] 


1  0  0 

0  cos(2^)  sin(2^) 

0  -sin(2^)  cos(2^) 


[V2(r)] 


cos(r) 
j  sin(r) 


y'sin(r) 

cos(r) 


[U3T(2t)\ 


cos(2r)  0  ysin(2r) 
0  1  0 
ysin(2r)  0  cos(2r) 


[v»] 


0 

+  iv 
e  J 


=>[U3T(2v)] 


cos(2v)  -y'sin(2v)  0 
-y  sin(2v)  cos(2v)  0 
0  0  1 


(29) 


(30) 


(31) 


This  special  unitary  elliptical  basis  changement  matrix  can  also  be  described  using  the  parameters  p  and  0 
which  correspond  to  the  polarisation  ratio  of  this  first  or  principal  Jones  vector  of  the  new  basis,  and  the 
matrix  [U3T\  is  thus  given  by  : 


[U3T(p,t)\ 


1 

Hd1 


cos(2£)  +  9t  e(p2e~2jt) 

-  j  sin(2£)  +  jZm(p2e~2M ) 
2j~sm(p) 


-  j sin(2£)  -  jZsnt(p2e  2,« ) 
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-  2  Me(p) 


2  jZ$m{pe  2j q 
291  e(pe-2J4) 


(32) 
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In  the  case  of  the  3x3  covariance  matrix  [C],  the  elliptical  polarisation  basis  changement  operator  is  given 
by: 


[ChPscfa.r.rftCpscfa.T.vft1  =  [U3Ap^)\c1u1c(p4)Y  (33) 


The  SU(3)  unitary  transformation  group  for  the  covariance  matrix  is  obtained  from  the  SU(3)  unitary 
transformation  group  for  the  coherency  matrix  with: 

[u3c  (2t,2T,2vj\  =  [u3  ][U3T  (2</>,2t,2v)\U3  ]"'  (34) 


where  : 


[v] 


i 

J2 


+  1  0  +1 
+  1  0  -1 
0  V2  0 


(35) 


Expressed  in  function  of  the  polarisation  ratio  parameters  (p  and  ^  )  the  special  unitary  elliptical  basis 
changement  matrix  [U3C\  can  also  be  written  following: 


[U,c(p)\ 


42pe-2Jt 

p2e2* 

■42  p 

H/f 

4~2p 

-42  pe+2ii 

e+2ii 

2.3  POLARIMETRIC  SAR  DATA  STATISTICS 


(36) 


A  polarimetric  radar  measures  the  complete  scattering  matrix  [S]  of  a  medium  at  a  given  incidence  angle  and 
for  a  given  frequency.  This  scattering  matrix  with  complex  elements  can  be  expressed  in  a  vector  form, 
following  : 


X  = 


4>S 


HV 


(37) 


where  the  factor  a I~2  arises  from  the  requirement  to  keep  the  norm  of  the  target  vector  X  an  invariant,  equal 
to  the  Frobenius  norm  (Span)  of  the  backscattering  matrix  [S],  namely  the  total  power  scattered  by  the  target. 
When  the  radar  illuminates  an  area  of  a  random  surface  of  many  elementary  scatterers,  the  one-look 
scattering  vector  X  can  be  modeled  as  having  a  multivariate  complex  Gaussian  probability  density  function 
[Lee  1999b],  with  : 


P{X) 


j_„-x'T[c]-'x 


w 


(38) 


where  [C]=E(XXT*)is  the  hermitian  covariance  matrix  of  the  scattering  vector  X. 

As  the  target  vector  k  is  constructed  from  a  linear  combination  of  the  elements  of  the  scattering  vector  X , 
following  : 


k  = 


1 

T2 


^HH  +  SVv 

i 

=  T2 

1 

0 

1 " 

Shh 

S  hh  ~  Syv 

1 

0 

-1 

42Shv 

_  2Shv 

_0 

2 

0  _ 

1 

Co 

_ 1 

(39) 


the  elements  of  the  target  vector  k  are  considered  as  having  the  same  complex  Gaussian  distribution  than  the 
elements  of  the  scattering  vector  X. 
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POLSAR  data  are  frequently  multi-look  processed  for  speckle  reduction,  or  data  compression.  The  relative 
polarimetric  information  is  thus  contained  in  the  expected  value  of  the  coherency  matrix  <[7]  >  representing 
the  spatial-averaged  distributed  target  given  by  : 


<m>  w 

iV  i=l  iy  i=l 

where  [Tt]  is  a  single-look  coherency  matrix  of  the  ith  pixel. 

It  has  been  shown  in  [Lee  1999b]  that  the  averaged  coherency  matrix  <[7]>  has  a  complex  Wishart 
distribution.  The  probability  density  function  is  given  by  : 


^(<[r]>/[r„]) 


^|([r])rVtr,(|rJ'<|r'>) 

n‘U^T{L)...T(L-p  +  l%TnX 


(41) 


where  [ Tm ]  is  global  coherency  matrix  with  [Tm\  =  E (k  kT*),  L  is  the  number  of  look  and  p  the  dimension  of 
the  target  vector  k ,  with p=3  for  the  reciprocal  case  (SHv=Svh)  and p=4  for  the  non-reciprocal  case. 

The  distribution  functions  for  dual  polarization  ( HH ,  VH),  (HV,  VV)  or  (HH,  VV)  can  be  derived  from  this 
complex  Wishart  distribution.  For  example,  if  only  complex  HH  and  VV  are  available,  p  =2,  and  for  single 
polarization,  p=l,  which  reduces  (9)  to  the  Chi-square  distribution  with  2 L  degree  of  freedom. 

For  the  dual  polarization  case  without  phase  difference  information  (\HH\,  \W\),  the  probability  density 

function  has  been  derived  [Lee  1994a].  Letting  Rj  =<  ||2  >  and  R2  =<  i^r  > ,  we  have: 


(n-l) 


nn+1  (RjR2 )  2  exp 


p(R,,R2)  = 


yi( Rj  /  Cjj  +  R2  /  C 22 ) 
l-\Pc  \2 


(C11C22)(n+I)/2Y(n)(l-\pc\2)\pc 


h-i 


2n 


R,R2  \pc 


C ijC 22  1  |  p 


c  I  J 


(42) 


where  In  (  )  is  the  modified  Bessel  function  of  the  nth  order,  Cn  =  and  C22  = 


3  SPECKLE  FILTERING 

3.1  Need  for  speckle  filtering 

Unlike  optical  remote  sensing  images,  characterized  by  very  neat  and  uniform  features,  SAR  images  are 
affected  by  speckle  [Goze  1993],  [Lee  1980],  [Lee  1981a],  [Lee  1981b],  [Lee  1983],  [Lee  1986a],  [Lee 
1986b],  [Lee  1994],  [Lee  1997],  [Lopes  1990],  [Lopes  1993],  [Touzi  1994].  Speckle  confers  to  SAR  images 
a  granular  aspect  with  random  spatial  variations.  Fig.  3  shows  an  example  of  single  polarization  speckled 
SAR  images. 
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Fig.  3:  Single-look  | .S', ,  | 2  (top)  and  |iS12 12  (bottom) images 


The  intensity  images  displayed  in  Fig.  3  show  a  poor  contrast,  as  well  as  a  random  aspect,  that  reduce  the 
possibilities  of  visual  interpretation  and  analysis  of  the  scene  under  consideration. 

The  discrimination  of  different  natural  media  by  comparing  intensity  to  a  fixed  threshold  leads,  in  general  to 
numerous  errors  due  to  the  high  variability  of  SAR  speckled  response.  Speckle  phenomenon  also  affects  the 
phase  of  scattering  coefficients  and  corrupts  polarimetric  information. 

The  image  of  Arg(Sn)  shown  in  Fig.  4  indicates  that  the  absolute  phase  of  a  scattering  coefficient  is  highly 
random  and  does  not  contain  evident  information.  Speckle  does  not  affect  similarly  different  polarimetric 
channels,  as  shown  in  the  between  channel  relative  phase  image,  Arg(SnS22)  and  in  the  color  coded 

image,  built  from  the  three  polarimetric  channel  intensities,  |*Sn|2,  \^n\2  and 

i?  =  |lS11|2,G  =  |1S12|2,5  =  |<S22|2 


Fig.  4:  Arg(Sll )  (top)  and  Arg{SuS22)  (middle)  images. 
Color  coded  image  R  =  |iSn|2?G  =  1^2 12, 5  =  |iS22|2  (bottom) 


Speckle  corrupts  polarimetric  observables  (phase  and  intensity)  in  an  important  way.  Specific  procedures 
have  to  be  used  to  retrieve  relevant  polarimetric  information  and  to  reduce  the  randomness  of  the  acquired 
signals. 
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3.2  Simple  speckle  model 


3.2.1  Single-polarization  multiplicative  speckle  model 

Speckle  confers  a  random  aspect  to  SAR  images,  but  may  not  be  considered  as  a  simple  noise.  It  is,  in  fact, 
tightly  related  to  SAR  measurement  principle.  [Goze  1993],  [Lee  1980],  [Lee  1981a],  [Lee  1981b],  [Lee 
1983],  [Lee  1986a],  [Lee  1986b],  [Lee  1994],  [Lee  1997],  [Lopes  1990],  [Lopes  1993],  [Touzi  1994]. 
Synthesized  SAR  data  may  be  considered  as  the  result  of  the  integration  of  a  scene  coherent  response  within 
each  resolution  cell,  resulting  from  the  convolution  of  the  SAR  impulse  response  with  the  coherent 
contribution  of  each  elementary  scatterer,  as  illustrated  in  Fig.  5.  As  the  number  of  contributing  scatterers, 
within  a  resolution  cell,  tends  to  be  large  (it  is  the  case  for  common  resolution  SAR  measurements),  the 
resulting  integrated  response  is  random  in  phase  and  amplitude  and  is  shown  to  follow,  over  homogeneous 
areas,  a  Normal  distribution. 


A  speckled  response  is  usually  represented  under  the  form  of  a  simple  product  model:  y  —  x.n  where 
y  represents  a  complex  speckled  scattering  coefficient,  x  the  original  unspeckled  scattering  coefficient  and 
n  the  multiplicative  speckle  contribution. 

The  speckle  term,  n  is  composed  of  independent  real  and  imaginary  parts,  following  both  real  centered 
Normal  distribution  Nc( 0, 1/2) .  The  corresponding  speckled  intensity,  7,  is 


7*  *  * 

=  yy  =xx  nn 


(43) 


=>  Y  =  X  nn 

Over  homogeneous  areas,  X  is  considered  to  be  constant  and  the  speckled  intensity  follows  an  exponential 
probability  density  function 


_Y 

p(Y)  =  X  V* 


(44) 


Its  first  two  moments  are  given  by: 


E(Y)  =  XE{nn)  =  X  with:  Var(Y)  =  X2  Var(nn)  =  X2  (45) 


3.2.2  POLARIMETRIC MULTIPLICATIVE  SPECKLE  MODEL 

This  speckle  model  may  be  extended  to  the  polarimetric  case  by  considering  that  polarimetric  channels  are 
affected  by  independent  multiplicative  speckle  components  : 
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<N  ' 
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Xn(nnnn) 

sa 
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*12 

nn 

= 

^12(^12^12) 

_^22_ 

*22  _ 

J1 22  _ 

^22  (^22^22) 

(46) 


One  may  note  that  the  multiplicative  assumption  is  in  general  not  valid  to  model  speckled  correlation  terms. 
3.3  Principle  of  scalar  speckle  filtering 


3.3.1  Incoherent  averaging 


As  presented  in  the  former  paragraph,  a  speckled  intensity,  7,  may  be  considered  as  a  random  variable 
whose  mean  value  equals  the  unspeckled  intensity,  X ,  but  affected  by  a  large  variance. 

The  principle  of  speckle  filtering  consist  in  reducing  the  variance  of  7  in  order  to  improve  the  estimate  of  its 
mean.  The  sample  mean,  (Y^ ,  is  defined  as  the  empirical  average  of  L  independent  realizations  of  a  speckled 
intensity  as  follows: 

(Y)=jiY,  <47> 

^  i=l 


It  can  be  shown  that,  over  homogeneous  areas,  this  estimate  of  X  follows  a  Gamma  density  function: 


P({Y)) 


llx{y)l-' 

T(L)Xl  e 


(48) 


and  has  the  following  two  first  moments 

E((Y))  =  X  with:  Var({Y))  =  X2  /  L  (49) 

It  is  possible  to  observe  from  (49)  that  as  the  number  of  independent  samples,  L,  reduces  to  1  the  variance  of 
the  estimate  intensity  increases,  whereas  incoherent  averaging  over  L  independent  realizations  permits  to 

reduce  the  variance  of  a  speckled  intensity  in  a  significant  way.  The  quantity  L  =  var(/)/i?(/)2  is  called  the 
Equivalent  Number  of  Looks  (ENL)  and  is  a  measure  of  speckle  importance. 

stdev(V) 


0  2  4  6  8  10 

rf&)) 

Fig.  6:  Occurrence  of  single-look  of  intensity  moments  over  a  homogeneous  areas 
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The  occurrence  plot  displayed  in  Fig.  6  clearly  shows  that  the  standard  deviation  and  the  mean  of  sampled 

intensities  are  linearly  related  over  homogeneous  areas.  The  slope  of  this  linear  relation  is  V Z  and  equals  1 
in  the  case  of  single-look  data  sets. 


3.3.2  Boxcar  filter 


The  boxcar  filter  is  a  direct  application  of  the  incoherent  averaging  described  by  (47)  to  the  case  of  an  image. 
Filtered  intensity  estimates,  X .  y. ,  are  constructed  by  computing  the  sample  mean  over  each  pixel 
neighborhood,  defined  by  a  sliding  window  of  (Nw  x  Nw )  pixels. 


X ,  =  (Z 


hj 


l’J  /  N, 


1 


Nw/2 


Nw/2 


i  m 


i+p,j+q 


p=-Nw/2q=-Nw/2 


(50) 


where  the  subscripts  i  and  j  correspond  to  the  considered  pixel  row  and  column  index  respectively.  Fig.  7 
shows  an  intensity  image  obtained  using  a  (7x7)  boxcar  filter.  This  images  shows  enhanced  contrast  and 
lower  random  aspect. 


Fig.  7:  |  filtered  image  using  a  boxcar  filter 


As  it  can  be  seen  in  Fig.  7,  the  boxcar  filter  is  characterized  by  two  main  limitations  : 

•  sharp  edges  are  generally  blurred 

•  point  scatterers  are  over  filtered  and  transformed  to  spread  targets 
Solutions  to  these  limitations  are  offered  by  the  refined  Lee  filter. 


3.3.3  J.S.  Lee  adaptive  filter 


J.  S.  Lee's  filter  determines  the  unspeckled  intensity  estimate  that  minimizes  a  mean  squared  error  [Lee 
1980],  [Lee  1981a],  [Lee  1981b],  [Lee  1983],  [Lee  1986a],  [Lee  1986b],  [Lee  1994]: 


( \x-x 


(51) 


This  MMSE  filter  is  based  on  a  linearized  speckle  model  leading  to  the  following  estimate  expression  : 

*=(yL.+*(My)J  <52> 

where  k  is  an  adaptive  filtering  coefficient,  based  on  local  statistics,  given  by  : 
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vn,(X)  var(r)-£!((r)K2 

var(r)  var<7}[l  +  <x„2]  '  ' 

2  1 

with  Gn  ——  the  a  priori  speckle  variance.  Over  homogeneous  areas,  var(X)  =  0  =>  k  =  0  and 
L 

X  =  (Y^n  ,  whereas  over  point  targets  and  highly  heterogeneous  areas,  k  =  1=>  X  =  Y  and  the  pixel 


intensity  remains  unaffected  by  the  filtering  procedure. 

In  order  to  reduce  the  sensitivity  of  the  adaptive  filtering  coefficient,  k ,  to  isolated  heterogeneities,  this  filter 
uses  directional  masks  to  determine  the  most  homogeneous  part  of  the  sliding  window  where  local  statistics 
have  to  be  estimated.  This  modification  permits  to  preserve  relatively  sharp  edges.  [Lee  1994],  [Lee  1997]. 


Fig.  8:  Examples  of  directional  masks 


■  0 

□  1 


The  Lee  filter  results  displayed  in  Fig.  9  demonstrate  the  effectiveness  of  this  adaptive  filtering  approach 


Fig.  9:  |  Sn^  filtered  image  using  a  (7x7)  Lee  filter 


3.4  Extension  to  the  polarimetric  case 

Speckle  filtering  is  based  on  incoherent  averaging  and  requires  to  handle  second  order  representations.  The 
intensity  information  used  in  the  scalar  case  has  to  be  extended  to  the  vector  case  when  dealing  with  two  or 
more  polarization  channels  in  order  to  take  into  account  the  different  intensities  as  well  as  the  cross¬ 
correlation  related  information.  A  simple  way  to  build  an  incoherent  polarimetric  representation  consists  in 
vectorizing  a  scattering  matrix  to  create  a  target  vector  and  computing  the  corresponding  (3x3)  covariance 
matrix,  [C3]  or  the  (3x3)  coherency  matrix,  [T3] . 


3. 4. 1  Polarimetric  boxcar  filter 

The  extension  of  the  boxcar  filter  to  the  polarimetric  case  is  straightforward.  The  estimated  covariance 
matrix  is  given  by  : 
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1 


Nw/2  N„,/2 


w  ^  w  p=-Nw/2q=-Nw/2 


J+q 


(54) 


where  the  subscripts  i  and  j  correspond  to  the  considered  pixel  row  and  column  index  respectively. 


3.4.2  J.  S.  Lee  polarimetric filter 

J.  S.  Lee  proposed  to  estimate  the  unspeckled  covariance  matrix  according  to  the  following  expression: 

[C]  =  ([C])+i([C]-([C]»  (55) 

where  k  remains  a  scalar  coefficient  computed  from  the  span  statistics,  span  =  Cn  +  C22  +  C33 . 

This  approximation  is  allows  to  filter  polarimetric  data  in  a  fast  and  simple  way  and  avoids  additional 
coupling  (or  cross-talk)  between  the  polarimetric  channels.  Fig.  10  shows  improved  color  coded  images 
processed  through  the  boxcar  and  J.  S.  Lee  filters. 


Fig.  10  :  Color-coded R  =  |iSn|2,G  =  |iS12|2,i?  =  |iS22 


filtered  image  using  a  (7x7)  boxcar  (top)  Lee's  filter  (bottom) 

3.5  Consequences  of  polarimetric  speckle  filtering 

It  was  seen  in  former  paragraphs  that  it  is  necessary  to  reduce  polarimetric  variables  random  aspect  by 
speckle  filtering  prior  to  any  interpretation  of  polarimetric  information.  The  incoherent  averaging  of  the 
coherency  [7]  or  covariance  [C]  matrices  has  an  important  impact  on  their  polarimetric  properties.  Speckle 
filtering  may  cause  a  loss  of  polarimetric  information  by  destroying  the  relation  between  and  [7]  or  [C] 
matrices. 

A  coherency  matrix  is  fully  defined  by  9  real  coefficients:  its  three  diagonal  terms  and  three  complex 
correlation  coefficients.  In  the  case  of  a  single  look  coherency  matrix,  all  three  correlation  coefficients  have 
unitary  modulus  and  one  of  their  phase  may  be  obtained  by  a  linear  combination  of  the  remaining  two, 
leaving  5  degrees  of  freedom.  A  relative  scattering  matrix  [Srel]  and  single-look  [7]  or  [C]  matrices  may  be 
related  in  a  unique  way  as  shown  in  the  following  example  : 
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(57) 


The  scattering  mechanism  may  then  be  interpreted  by  comparing  [Srel  ]  to  canonical  examples.  After  speckle 

filtering,  this  may  not  be  true  anymore.  In  a  general  case,  the  modulus  of  correlation  coefficients  is  inferior  to 
one  and  the  phase  terms  are  linearly  independent. 


SVSU 


<(\s„ 


& 


and 


(58) 


In  such  a  case,  the  coherency  matrix  is  said  to  be  distributed  and  cannot  be  related  to  a  coherent  scattering 
matrix. 


Fig.  11:  Argument  (top)  and  modulus  (bottom)  of 


after  application  of  a  Lee  filter 


The  correlation  coefficient  displayed  in  Fig.  1 1  shows  a  varying  modulus  over  the  selected  scene,  indicating 
that  the  degree  of  correlation  might  be  related  to  the  nature  of  the  scattering  medium.  The  additional 
information  contained  in  the  cross-correlation  terms  will  be  exploited  by  incoherent  decomposition 
theorems  to  extract  even  more  characteristics  from  polarimetric  data  sets. 


4  POLARIMETRIC  TARGET  DECOMPOSITION  THEOREMS 

There  is  currently  a  great  deal  of  interest  in  the  use  of  polarimetry  for  radar  remote  sensing.  In  this  context, 
an  important  objective  is  to  extract  physical  information  from  the  observed  scattering  of  microwaves  by 
surface  and  volume  structures.  The  most  important  observable  measured  by  such  radar  systems  is  the  3x3 
coherency  matrix  [7].  This  matrix  accounts  for  local  variations  in  the  scattering  matrix  and  is  the  lowest 
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order  operator  suitable  to  extract  polarimetric  parameters  for  distributed  scatterers  in  the  presence  of  additive 
(system)  and/or  multiplicative  (speckle)  noise. 

Many  targets  of  interest  in  radar  remote  sensing  require  a  multivariate  statistical  description  due  to  the 
combination  of  coherent  speckle  noise  and  random  vector  scattering  effects  from  surface  and  volume.  For 
such  targets,  it  is  of  interest  to  generate  the  concept  of  an  average  or  dominant  scattering  mechanism  for  the 
purposes  of  classification  or  inversion  of  scattering  data.  This  averaging  process  leads  to  the  concept  of  the 
«  distributed  target »  which  has  its  own  structure,  in  opposition  to  the  stationary  target  or  « pure  single 
target »  [Huynen  1970][Pottier  1992]. 

Target  Decomposition  theorems  are  aimed  at  providing  such  an  interpretation  based  on  sensible  physical 
constraints  such  as  the  average  target  being  invariant  to  changes  in  wave  polarization  basis. 

Target  Decomposition  theorems  were  first  formalized  by  J.R.  Huynen  but  have  their  roots  in  the  work  of 
Chandrasekhar  on  light  scattering  by  small  anisotropic  particles.  Since  this  original  work,  there  have  been 
many  other  proposed  decompositions.  We  classify  four  main  types  of  theorem: 

1 .  Those  employing  coherent  decomposition  of  the  scattering  matrix  (Krogager,  Cameron). 

2.  Those  based  on  the  dichotomy  of  the  Kennaugh  matrix  (Huynen,  Barnes). 

3.  Those  based  on  a  “model-based”  decomposition  of  the  covariance  or  the  coherency  matrix 
(Freeman  and  Durden,  Dong). 

4.  Those  using  an  eigenvector  /  eigenvalues  analysis  of  the  covariance  or  coherency  matrix 
(Cloude,  VanZyl,  Cloude  and  Pottier). 

A  complete  description  of  all  these  different  Polarimetric  Target  Decomposition  can  be  found  in  [Cloude 
1996],  and  we  focus  here  on  the  H  /  A  /  a  decomposition  theorem  which  will  be  used  further  in  the 
polarimetric  classifications. 

4.1  The  H  /  A  /  a  Polarimetric  Decomposition  Theorem 

In  1997,  S.R.  Cloude  and  E.  Pottier  proposed  a  method  for  extracting  average  parameters  from  experimental 
data  using  a  smoothing  algorithm  based  on  second  order  statistics  [Cloude  1997].  This  method  does  not  rely 
on  the  assumption  of  a  particular  underlying  statistical  distribution  and  so  is  free  of  the  physical  constraints 
imposed  by  such  multivariate  models.  An  eigenvector  analysis  of  the  3x3  coherency  matrix  [7]  is  used  since 
it  provides  a  basis  invariant  description  of  the  scatterer  and  also  provides  a  decomposition  into  types  of 
scattering  process  (the  eigenvectors)  and  their  relative  magnitudes  (the  eigenvalues).  This  original  method, 
based  on  an  eigenvalue  analysis  of  the  coherency  matrix,  employs  a  3 -level  Bernoulli  statistical  model  to 
generate  estimates  of  the  average  target  scattering  matrix  parameters.  This  alternative  statistical  model  sets 
out  with  the  assumption  that  there  is  always  a  dominant  'average'  scattering  mechanism  in  each  cell  and  then 
undertakes  the  task  of  finding  the  parameters  of  this  average  component  [Cloude  1997]. 


4. 1. 1  Eigenvector-based  Decomposition. 

The  instantaneous  (single  pixel)  target  return  from  a  spatially  extended  target  can  be  characterized  either  by 
its  complex  scattering  matrix  [S]  which  relates  to  received  spatial-voltage,  or  by  its  3x3  coherency  matrix  [7] 
which  relates  to  spatial-power.  In  the  case  of  spatial-averaging,  it  is  customary  to  consider  the  expected  value 
of  the  coherency  matrix  <[7]>  as  representing  the  averaged  distributed  target,  as  : 

From  this  estimate,  the  eigenvectors  and  eigenvalues  of  the  3x3  hermitian  coherency  matrix  <[7]>  can  be 
calculated  to  generate  a  diagonal  form  of  the  coherency  matrix  which  can  be  physically  interpreted  as 
statistical  independence  between  a  set  of  target  vectors  [Cloude  1992][  Cloude  1996].  The  coherency  matrix 
<[7]>  can  be  written  in  the  form  of: 
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([T}=[u3Mu3r 


(60) 


where  [i7]  is  a  3x3  diagonal  matrix  with  nonnegative  real  elements,  and  [t/3]  =  \uj  u2  ]  is  a  3x3 

unitary  matrix  of  the  SU(3)  group,  where  uu  M2,  and  u3  are  the  three  unit  orthogonal  eigenvectors. 

By  finding  the  eigenvectors  of  the  3x3  hermitian  coherency  matrix  <[7]>,  such  a  set  of  3  uncorrelated  targets 
can  be  obtained  and  hence  a  simple  statistical  model  can  be  constructed,  consisting  of  the  expansion  of  <[7]> 
into  the  sum  of  3  independent  targets,  each  of  which,  represented  by  a  single  scattering  matrix.  This 
decomposition  can  be  written  following: 


(M)=Z47;]=IA",-K:r  («D 

i=l  i=l 

where  the  real  numbers  are  the  eigenvalues  of  <[7]>  and  represent  statistical  weights  for  the  three 
normalized  component  targets  [7}]  [Cloude  1992]. 

If  only  one  eigenvalue  is  nonzero  then  the  coherency  matrix  <[7]>  corresponds  to  a  pure  target  and  can  be 
related  to  a  single  scattering  matrix.  On  the  other  hand,  if  all  eigenvalues  are  equal,  the  coherency  matrix 
<[7]>  is  composed  of  three  orthogonal  scattering  mechanisms  with  equal  amplitudes,  the  target  is  said 
random  and  there  is  no  correlated  polarized  structure  at  all. 

Between  these  two  extrema,  there  exists  the  case  of  partial  targets  and  where  the  coherency  matrix  <[7]>  has 
non-zero  and  non-equal  eigenvalues.  The  analysis  of  its  polarimetric  properties  requires  a  study  of  the 
eigenvalues  distribution  as  well  as  a  characterization  of  each  scattering  mechanism  of  the  expansion. 

To  introduce  the  degree  of  statistical  disorder  of  each  target,  the  entropy  (77)  is  defined  in  the  Von  Neumann 
sense  from  the  logarithmic  sum  of  eigenvalues  of  <[7]>  [Cloude  1996]  [Cloude  1997],  as: 

H  =  -YJPllog3(Pi)  (62) 

i=l 


where  PL  are  the  probabilities  obtained  from  the  eigenvalues  X[  of  <[7]>  with: 


P  = 


1 


1  j=3 

2a 


j=l 


(63) 


If  the  entropy  H  is  low  then  the  system  may  be  considered  as  weakly  depolarizing  and  the  dominant  target 
scattering  matrix  component  can  be  extracted  as  the  eigenvector  corresponding  to  the  largest  eigenvalue  and 
ignore  the  other  eigenvector  components. 

If  the  entropy  H  is  high  then  the  target  is  depolarizing  and  we  can  no  longer  consider  it  as  having  a  single 
equivalent  scattering  matrix.  The  full  eigenvalue  spectrum  must  be  considered. 

Further,  as  the  entropy  H  increases,  the  number  of  distinguishable  classes  identifiable  from  polarimetric 
observations  is  reduced.  In  the  limit  case,  when  H=1 ,  the  polarization  information  becomes  zero  and  the 
target  scattering  is  truly  a  random  noise  process  . 


While  the  entropy  is  a  useful  scalar  descriptor  of  the  randomness  of  the  scattering  problem,  it  is  not  a  unique 
function  of  the  eigenvalue  ratios.  Hence,  another  eigenvalue  parameter  defined  as  the  anisotropy  A  can  be 
introduced,  with  : 


~ 

X2  +  A3 


(64) 
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When  A=0  the  second  and  third  eigenvalues  are  equal.  The  anisotropy  may  reach  such  a  value  for  a 
dominant  scattering  mechanism,  where  the  second  and  third  eigenvalues  are  close  to  zero,  or  for  the  case  of  a 
random  scattering  type  where  the  three  eigenvalues  are  equal. 

The  condition  for  <[7]>  to  have  such  an  equivalent  scattering  matrix  [S]  is  for  both  the  target  entropy  H  and 
the  anisotropy  A  to  be  equal  to  zero,  which  corresponds  to  a  single  nonzero  eigenvalue  (k\)  [Cloude 
1992][Cloude  1996].  In  this  case  the  coherency  matrix  <[7]>  has  rank  r=l ,  and  can  be  expressed  as  the  outer 
product  of  a  single  target  vector  k /  with: 

([r])  =  • k]T  =  A jUj  -u]T  (65) 


where  \  —  2yA0  +  B0J  is  equal  to  the  Frobenius  norm  (Span)  of  the  corresponding  scattering  matrix,  and 
where  the  corresponding  unit  target  vector  is  expressed  as  follows  : 


u/ 


Aa0x, 


2A0 
C  +  jD 
H-jG 


jB~^Be+jarctan(D/C) 
_  g  ^-j  arctan( G / H) 


(66) 


It  is  interesting  to  note  that  the  modulus  of  the  three  components  of  the  unit  target  vector  are  directly  function 
of  the  three  «  Huynen  target  generators  ». 

Without  using  ground  truth  measurements,  this  polarimetric  parameterisation  of  the  unit  target  vector  u 
involves  the  fit  of  a  combination  of  three  simple  scattering  mechanisms  :  surface  scattering,  dihedral 
scattering  and  volume  scattering,  which  are  characterized  from  the  three  components  (target  generators)  of 
the  unit  target  vector  such  as: 

Surface  scattering:  A0  »  B0+B  ,  B0  -  B 

Dihedral  scattering:  B0  +  B  »  A0  ,  B0  -  B 

Volume  scattering:  B0  -  B  »  A0  ,  B0+B 


4.1.2  Probabilistic  model  for  random  media  scattering. 

In  previous  publications  [Cloude  1995]  [Cloude  1996]  [Cloude  1997],  a  parameterisation  of  the  eigenvectors 
of  the  3x3  coherency  matrix  [7]  has  been  introduced  for  the  case  of  scattering  medium  which  does  not  have 
azimuthal  symmetry  [Nghiem  1992],  and  which  takes  the  form  shown  in  (67). 


u  = 


cos  a  sin  a  cos  / 3  sin  a  sin  (5 


(67) 


It  follows  a  revised  parameterisation  of  the  coherency  matrix  [Cloude  1997],  as 


([r])  =  [t/j] 


Aj  0 

0  A , 


0 

0 


0  0  A, 


M' 


(68) 
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sina3cos(33e 

sin  a,  sin 

sin  a2  sin  fi2e^2  + 

^  sin  a3  sin  (53eJ 

(<%  +  <4) 
foi  +  k) 


(69) 


The  parameterisation  of  a  3x3  unitary  matrix  [U3\  in  terms  of  column  vectors  with  different  parameters  a1  P1 
etc...  is  made  so  as  to  enable  a  probabilistic  interpretation  of  the  scattering  process.  In  general,  the  columns 
of  the  3x3  unitary  matrix  are  not  only  unitary  but  mutually  orthogonal.  This  means  that  in  practice  a1  a2  and 
a3  are  not  independent. 

In  this  case  a  statistical  model  of  the  scatterer  is  considered  as  a  3  symbol  Bernoulli  process  i.e.  the  target  is 
modeled  as  the  sum  of  three  [S]  matrices,  represented  by  the  columns  of  [ U3 ]  in  (69),  occuring  with 
probabilities  P. , given  from  (15)  by  the  normalized  eigenvalues  so  that  Pj  +  P2  +  P3  =  1  [Cloude  1997]. 

In  this  way  for  example,  the  parameter  a  follows  a  random  sequence: 


a 


-{ 


a,  a 2  a 3  a 2  a 3  a 2  a 3a,a 


'/•••} 


(70) 


and  the  best  estimate  of  this  parameter  is  given  by  the  mean  of  this  sequence,  easily  evaluated  as  [Cloude 
1997] : 

a  =  PjCC,  +  P2a2  +  P3a3  (71) 


In  this  way,  the  mean  parameters  of  the  dominant  scattering  mechanism  are  extracted  from  the  3x3 
coherency  matrix  as  a  mean  target  vector  u0,  such  that: 


u0=e 


j<t> 


cosa 


JS 


sinacosfie 


.  R  Jr 

sin  a  sin  pe  — 


(72) 


where  the  parameters  a  j3  8  and  jare  defined  in  (19),  and  where  $  is  physically  equivalent  to  an  absolute 
phase. 


4 . 1. 3  The  roll  invarlance  property. 


One  of  the  most  important  property  in  Radar  Polarimetry  concerns  the  roll  invariance.  The  effect  of  rotation 
around  the  radar  line  of  sight  [Cloude  1996]  can  be  generated  as: 

([J(0)])  =  [u/]([r])[t/;]"7  (73) 

where  I U3  I  is  the  unitary  similarity  rotation  matrix,  given  by: 


1  0  0 
0  cos  20  sin  2d 
0  -sin  26  cos 26 


(74) 
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According  to  the  eigenvector-based  decomposition  approach,  the  coherency  matrix  can  be  written  in  the 
form  of: 


*  T 


(75) 


COSOCj 

cosa2eJ 

j(f)'3 

cosa3eJ  Ti 

sin  a,cosP'e^' 

sina2cosfi'e^2  +  ^ 

sinatCOSplM*^ 

i  y f 

sin  a )  sin  f3I'eJ  7 1 

sin  a2  sin  fi2e^2  +  - 

1  since,  sin  fiAr’>+ $) 

(76) 


where  [2]  is  the  same  3x3  diagonal  matrix  with  nonnegative  real  elements,  [t/j]  =  [y7  v2  v3  is  the 

new  3x3  unitary  matrix  of  the  SU(3)  group,  where  yi,  y2,  and  v3  are  the  new  three  unit  orthogonal 
eigenvectors. 

Following  the  parameterisation  of  the  3x3  unitary  matrix  \U  3],  it  can  be  seen  that  the  three  parameters  a2 
and  a3  remain  invariant,  as  the  three  eigenvalues  (Aj  A2  A3). 


It  follows  that  the  mean  parameter  a ,  and  the  two  important  scalar  functions  of  the  eigenvalues,  the  entropy 
H  and  the  anisotropy  A,  are  three  roll-invariant  parameters. 

Among  the  mean  parameters  (a,  S  and  y)  of  the  dominant  scattering  mechanism  which  can  be  extracted 
from  the  3x3  coherency  matrix,  it  is  now  clear  from  the  above  analysis,  that  for  random  media  problems,  the 
main  parameter  for  identifying  the  dominant  scattering  mechanism  is  a.  The  three  others  parameters  (/£  8 
and  y)  can  be  used  to  define  the  target  polarisation  orientation  angle  [Pottier  1998][Pottier  1999] [Schuler 
1999]. 

In  previous  publication  [Cloude  1997],  it  has  been  shown  that  the  useful  range  of  the  parameter  a 
corresponds  to  a  continuous  change  from  surface  scattering  in  the  geometrical  optics  limit  ( a=0° )  through 
surface  scattering  under  physical  optics  to  the  Bragg  surface  model,  encompassing  dipole  scattering  or  single 
scattering  by  a  cloud  of  anisotropic  particles  ( a=45° ),  moving  into  double  bounce  scattering  mechanisms 
between  two  dielectric  surfaces  and  finally  reaching  dihedral  scatter  from  metallic  surfaces  ( a=90° ). 

The  a  parameter  estimate  is  related  directly  to  underlying  average  physical  scattering  mechanism,  and  hence 
may  be  used  to  associate  observables  with  physical  properties  of  the  medium  [Cloude  1997].  Fig.  12  shows 
these  three  roll-invariant  parameters  :  H ,  A  and  a. 


0  0.5  1  0  0.5  1  0°  45°  90° 

Entropy  H  Anisotropy  A  a  parameter 

Fig.  12  :  Roll-invariant  parameters  :  H,  A  and  a 
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The  figure  Fig.  12c  corresponding  to  the  a  parameter  image  shows  that  low  value  occurs  over  the  ocean 
region,  indicative  of  dominant  single  scattering  (a  =  0°).  Urban  area  and  parkland  areas  consist  of  medium 
and  high  a  parameter  values  (45 °<  a  <90°). 

The  figure  Fig.  12a  corresponding  to  the  entropy  image  shows  that  low  entropy  scattering  occurs  over  the 
ocean  (scattering  by  a  slightly  rough  surface).  High  entropy  occurs  over  the  parkland  areas.  At  this  resolution 
,  the  urban  area  consists  of  a  mixture  of  low  and  high  entropy  processes,  which  are  due  to  the  different 
street/building  classes  which  are  aligned  with  the  radar  look  direction,  or  aligned  somewhat  off  boresight,  or 
45°  aligned. 

The  figure  Fig.  12b  corresponding  to  the  anisotropy  image  shows  that  low  anisotropy  scattering  occurs  both 
over  the  ocean  region  and  parkland  areas.  The  fact  that  the  second  and  third  eigenvalues  are  equal, 
corresponds  either  to  a  single  dominant  scattering  mechanism  or  to  a  random  scattering  type.  The  urban  area 
and  the  coastal  sea  consist  of  a  mixture  of  medium  and  high  anisotropy  (presence  of  a  second  mechanism). 
Although  we  can  already  identify  several  clear  classes  based  on  entropy  alone,  further  information  can  be 
gleaned  from  the  angle  a  and  from  the  anisotropy  A,  which  can  distinguish  between  the  high  entropy  park 
and  urban  environments  since  the  latter  contains  (moderate  to  high  entropy)  dihedral  scattering. 


4.1.4  The  three-dimensional  H/A  /  a  classification  space. 

In  previous  publication  [Cloude  1997],  an  unsupervised  classification  scheme  has  been  introduced,  based  on 
the  use  of  the  two-dimensional  H  /  a  classification  plane,  where  all  random  scattering  mechanisms  can  be 
represented.  The  key  idea  is  that  entropy  arises  as  a  natural  measure  of  the  inherent  reversibility  of  the 
scattering  data  and  that  the  alpha  angle  (a)  can  be  used  to  identify  the  underlying  average  scattering 
mechanisms. 

The  H  /  a  classification  plane  is  sub-divided  into  nine  basic  zones  characteristic  of  classes  of  different 
scattering  behavior,  in  order  to  separate  the  data  into  basic  scattering  mechanisms,  as  shown  on  Fig.  13a.  The 
location  of  the  boundaries  within  the  feasible  combinations  of  H  and  a  values  is  arbitrary  and  generically, 
i.e.  based  on  the  general  properties  of  the  scattering  mechanisms.  There  is  of  course  some  degree  of 
arbitrariness  on  the  setting  of  these  boundaries  which  are  not  dependent  on  a  particular  data  set.  This 
segmentation  of  the  H  /  a  classification  plane  is  offered  merely  to  illustrate  the  unsupervised  classification 
strategy  and  to  emphasize  the  geometrical  segmentation  of  physical  scattering  processes. 

Complete  details  of  the  physical  scattering  characteristics  of  each  of  the  nine  zones  can  be  found  in  [Cloude 
1997][Pottier  2000]. 

The  distribution  of  the  San  Francisco  bay  POLSAR  data  on  the  H  /  a  classification  plane  is  shown  on  Fig. 
13b. 


Low 


Medium  Entropy  High  Entropy 


Fig.  13a  :  The  H  /  a  classification  plane. 
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Fig.  13b  :  Polsar  data  distribution  in  the  H  /  a 
classification  plane. 


RTO-EN-SET-081 


4-21 


Advanced  Concepts  in  Polarimetry  -  Part  1 


ORGANIZATION 


Inherent  of  the  spatial  averaging,  the  entropy  H  may  increase,  and  the  number  of  distinguishable  classes 
identifiable  from  polarimetric  observations  is  reduced.  For  example,  the  feasible  region  of  the  H  /  a 
classification  plane  is  rapidly  shrinking  for  high  values  of  entropy  (H=0.9),  where  ^  parameter  reaches  the 
limited  value  of  60°. 

A  high  entropy,  H=0.9 ,  may  correspond  to  two  limit  types  of  scattering  process  with  associated  eigenvalues 
spectra  given,  for  example,  by  :  A\=l,  A2 =0.4,  A3=0.4  or  by  :  A\  =1,  A2=l,  M =0.3. 

To  distinguish  between  these  two  different  types  of  scattering  process,  it  is  necessary  to  introduce  the 
anisotropy  information,  which  takes  the  corresponding  values  A=0  and  A=0.54  for  the  two  previous 
examples. 

In  order  to  extend  the  classification  scheme  and  to  improve  the  capability  to  distinguish  different  types  of 
scattering  process,  it  is  proposed  to  use  some  combinations  between  entropy  and  anisotropy  information,  as 
shown  on  Fig.  14.  The  .*  operation  represents  the  element  by  element  multiplication  of  two  matrices. 

The  examination  of  the  different  figures  corresponding  to  the  different  combinations  between  entropy  and 
anisotropy  images  leads  to  the  following  remarks  : 

1)  -  The  (1-H)(1-A)  image  corresponds  to  the  presence  of  a  single  dominant  scattering  process  (low 
entropy  and  low  anisotropy  with  A2  ~  A 3  ~0). 

2)  -  The  H(l-A)  image  characterizes  a  random  scattering  process  (high  entropy  and  low  anisotropy 
with  A2  ~  ^3  &AJ 1) 

3)  -  The  HA  image  relates  to  the  presence  of  two  scattering  mechanisms  with  the  same  probability 
(high  entropy  and  high  anisotropy  with  A3  &0). 

4)  -  The  (1  -H)A  image  corresponds  to  the  presence  of  two  scattering  mechanisms  with  a  dominant 
process  (low  to  medium  entropy)  and  a  second  one  with  medium  probability  (high  anisotropy  with  A3  &0). 


These  remarks  are  confirmed  by  the  analysis  of  the  distribution  of  the  San-Francisco  bay  POLSAR  data  in 
the  extended  and  complemented  three-dimensional  H  /  A  /  a  classification  space,  as  shown  on  Fig.  6  This 
representation  shows  that  it  is  possible  to  discriminate  new  classes  using  the  anisotropy  value. 

For  example,  it  is  now  possible  to  notice  that  there  exists  in  the  «  Low  Entropy  Surface  Scattering  »  area 
(Z9)  a  second  class  associated  with  a  high  anisotropy  value  and  which  corresponds  to  the  presence  of  a 
second  physical  mechanism  which  is  not  negligible. 

Identical  remarks  can  be  made  concerning  the  «  Medium  Entropy  Vegetation  Scattering  »  area  (Z5)  and  the 
«  Medium  Entropy  Multiple  Scattering  »  area  (Z4).  Due  to  the  spread  of  the  POLSAR  data  along  the 
anisotropy  axis,  it  is  now  possible  to  improve  the  capability  to  distinguish  different  types  of  scattering 
process  which  have  quite  the  same  high  entropy  value: 

-  High  entropy  and  low  anisotropy  correspond  to  random  scattering. 

-  High  entropy  and  high  anisotropy  correspond  to  the  presence  of  two  scattering  mechanisms  with 
the  same  probability. 

From  the  analysis  of  the  different  images  shown  on  Fig.  14  and  from  the  distribution  of  the  San-Francisco 
bay  POLSAR  data  in  the  H  /  A  /  a  classification  space  shown  on  Fig.  15.,  we  can  conclude  that  the 
anisotropy  has  to  be  considered  now  as  a  key  parameter  in  the  polarimetric  analysis  and/or  inversion  of 
POLSAR  data. 

The  information  contained  in  these  three  roll-invariant  parameters  extracted  from  the  local  estimate  of  the 
3x3  hermitian  coherency  matrix  <[7]>,  corresponds  to  the  type  of  scattering  process  which  occurs  within  the 
pixel  to  be  classified  (combination  of  entropy  H  and  anisotropy  A)  and  to  the  corresponding  physical 
scattering  mechanism  (^  parameter). 
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Fig.  14  :  Combinations  between  entropy  and  anisotropy  images. 
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Fig.  15  :  Distribution  of  the  San-Francisco  bay  POLSAR  data  in  the 
H  /  A  /  a  classification  space 
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